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In terms of effective field theory and mixed propagator approach, we show that there is a larger 
hidden effect of isospin breaking in p° — > 7r°7 decay due to an ui exchange, p° — > lu — *■ Tr 0r y. The 
branching ratio is predicted as B(p° — * ir°'y) = (11.67±2.0) x 10 -4 , which is much larger than Partcile 
Data Group's datum (6.8±1.7) x 10~ 4 and one of charged mode, B(p ± -> tt ± 7) = (4.5±0.5) x 10 -4 . 
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I. INTRODUCTION 

In hadron physics, the isospin symmetry or charge symmetryjlj is broken by inequality of the light quark masses, 
especially m u ^ m,d, and electromagnetic interaction of hadrons. This breaking of isospin symmetry induces various 
measurable physics effects such as ir° — r\, A — S° mixing, oj — » 7r + 7r~ decay, etc. It is common knowledge that 
the effects of isospin breaking can be omittd in the isospin conservation processes. However, this argument is not 
always right. The purpose of this letter is to show that a large isospin breaking effect is in isospin conservation decay 
p° — > 7r°7 (so-called hidden isospin symmetry breaking effect). Thus the real branching ratio for this decay should be 
much larger than datum cited by PDG-2000Q and charged mode. 

The anomalous-like radiative decays of light flavour vector mesons have been observed by several group: The 
branching ratio for charged p is B(p ± — > 7r ± 7) = (4.5 ± 0.5) x 10~ 4 (3|. The branching ratios for p° and uj decays 
were first obtained in refs.Jf], |5[ || by using the data of e + e~ collider in Neutral Detector, but the triangle anomaly 
contribution of QCD was ignored in their analysis. Benayoun et al. has reanalyzed the p° and u> decay to a pseudoscalar 
meson plus a photon via taking into account the triangle anomaly contribution 0. The authors use two models (Mi 
and M2)||, |j| to fit data. It is a bit surprising that they found two local minima for both two model fits, one with 
X 2 /Ndf — 0.5 (solution A) and another one with \ 2 l^df — 0-7 (solution B). To see table 1, where phase <!>y of vector 
mesons is defined via e + e~ — * tt°j cross section (equation (1) of ref.Q). Although the solution A has smaller x 2 /Ndf, 
and the phase difference $^ — <f> p is found around 210° in the solution A, close to expectation from the quark model 
(180°), the p° branching ratio is much larger than expected from the charged mode or from SU(3) symmetry. Thus the 
authors took the solution B as final result for agreeing with isospin symmetry. However, this conclusion is incorrect. 
In this paper let us see what happens there. 

The paper is organized as follows: In section 2, we first give the formula of the transition amplitude for p, oj — > ttj. 
Then we will prove this formula via two independent methods: effective field theory approach and mixed propagator 
approach. In section 3, we will provide final numeric results, and a brief conclusion. 
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TABLE I: Branching ratios and phases obtained from fit to e + e — > 7r°7 cross section using two p° models and assuming 
$p = $w The table is abstracted from table 4 of ref.[7]. 

r(7T° -» 77 )=7.7eV Model Mi Model M 2 



sol. A sol. B sol. A sol. B 



p° -» 7r°7 (units 10" 4 ) (11.51 ± 2.00) (6.17 ± 1.57) (11.67 ± 2.00) (6.77 ± 1.72) 

cu -> 7r°7 (units 10" 2 ) (8.62 ± 0.26) (8.39 ± 0.24) (8.64 ± 0.27) (8.39 ± 0.24) 

-> 7T°7 (units 10~ 3 ) (1.15 ±0.13) (1.20 ±0.16) (1.21 ±0.16) (1.26 ±0.17) 

$ P = $ w [degrees] -94 ± 8 124 ± 9 -90 ± 7 125 ±11 

$ [degrees] 114 ± 14 248 ± 9 119+^ 248 + ^ 

X 2 / dof 29/57 38/57 29/57 36/57 



II. TRANSITION AMPLITUDE FOR V -> tt 7 

The transition matrix element for vector mesons decay to a pion and a photon is 

< 7T(k) 7 ( qi )\V(q 2 ) >= ifv.je^q^e^ep, (1) 

where e M and are polarized vector of vector mesons and photon respectively. At the large N c limit and chiral limit, 
if we ignore possible contribution from resonance exchange, the exact isospin symmetry implies 

f± - f (0) -if (0) (2) 

Jp ± TT ± f — J p 7r 7 — r^J^ 7T 1 , 

where the superscript "(0)" denotes absence of the resonance exchange. A complete consideration for p° and u>° 
decays should include the contribution from the following resonance exchange, 

( c ) ,(o) , n p ^(mp)/^ 7r ) o 7 



V ' = f 

• / p°7r°7 • / p°ir°7 ' ™2 

""p 

Ac) = AO) n p^(^)/pQ OT o 7 

Juno-y J ^°-y^ ml-m 2 p + im p T p ' { > 

where the momentum-dependent p° — uj mixing amplitude H puJ (q 2 ) is defined in the following effective lagrangian. 
The relations in equation (||) holds more generally. In the other words, the couplings /jo^.o 7 and ff® can include 
corrections beyond large N c limit and chiral limit. Meanwhile, the relation ^ will be broken, but this breaking is 
slight. This point can be checked by using datum of ( o ± — * 7r ± 7 decay and to — > 7r°7 decay if we believe / p ± 7r ± 7 = f ^\ 0l - 
It must be point out that the equation (^]) is a non-perurbative results instead of a perturbative expression. In the 
rest of this section we will prove it by two independent methods: effective field theory approach and mixed propagator 
approach, respectively. 



A. Effective field theory approach 

The most general effective lagrangian concerning p° — > 7r°7 and to — > 7r°7 decays as follow 

£ = £q + £l + £2 + -CchPT + £wzw, (4) 

where £chPT is lagrangian of chiral perturbative theory, £wzw is Wess-Zumino-Witten lagrangian, £0 and £1 arc 
lagrangians of free fields and interaction for vector mesons respectively, £2 is counterterm lagrangian for £q and £1 
Em . Explicitly, £q and £1 can be written as follows 
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^0 = - 1 (dp,pi-d uP i)(d»p i »-d»p i »)-~(d^ u -d^)(dw-d»un 



~ 2 i in 1 ~ 2 u 
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+e^e«V W 3 7 r(g 2 K(?)M^ 

+e^f bJ ^(q 2 )q^(q)d a A fJ (x)TT°(x) + •••}. (5) 

Several remarks are related to lagrangian (^): 

1. Focusing on any low energy effective theory, a general knowledge is that any coupling should be momentum- 
dependent (TT). In lagrangian (||) this truth has been exhibited via momentum-dependence of form factors 
fpinr(q 2 ), fui3ir(q 2 ), etc.. All these form factors should be real function of vector meson four-momentum square, 
q 2 . So that the unitarity of S'-matrix prevents us from taking complex mass square in on-shell transition 
amplitude. 

2. rhp and rh u in Cq are not physical masses of p° and uj, since they will be shifted due to p° — uj mixing jnj]. The 
physical masses of p° and uj correspond to poles in their type I complete propagators (fig. 1) [ po| , in which we 
have considered contribution from p° — uj mixing completely: 

a!1„(<7 2 ) - 



A (I) (J} = ~ ig ^ (Q) 



= H £ s + & 8 s s + ... 

— _| ^ x _|_ ^ g _g x + . . . 

FIG. 1: Non-perturbative correction of p° — uj mixing to propagators of p° and uj. The 
physical masses of p° and uj correspond to the poles of type I complete propagators. Here 
solid lines and dash lines denote nonphysical propagator of p° and uj respectively, and bold 
solid line and bold dash line denote type I complete propagators of p° and uj respectively. 

These type I complete propagators give the physical masses of p° and u> as follows 

m 2 p = m 2 + H pu (m 2 ), ml = ml - U pUJ (ml). (7) 

It is well-known that physical mass obtained from the pole of complete propagator should be equal to one 
obtained from orthogonal rotating p° and uj to their mass eigenstates. 

3. When we write interaction lagrangian ([&]), we have conveniently assumed that widths of p and uj are generated 
dynamically by pion loops. This assumption can be released, and it does not affect the following formal 
discussion. 

In Heisenberg picture, the transition matrix ([!]) can be expressed in terms of LSZ reducation formula [|To| 

< n(k)j(qi),out\V(q 2 ),in > 

x(d 2 x + m l)d 2 (d 2 + m 2 p )(0\T{n(x)A,(y)VAz)}\Oh, (8) 

where Z3 are renormalization constant of wave function, and the subscript "H" denotes in Heisenberg picture. The 
above expression can be transformed to interaction picture via 

(0|T{7r(x)^(y)K(2)}|0) H =< 0|T{ 7 r( 2 ;)A AI (y)V;(^e l / !i4a; ' (£7+£chPT+£wzw)(a; ' ) }|0 > : , (9) 

where the subscript "I" denotes in interaction picture. 

In the language of Feymann diagram, the non-perturbative results can be obtained by summing all diagrams of 
perturbative expansion (fig. 2). Here we focus on p° — > 7r°7 decay, and the discussion for to — > ir°j is similar. In fig. 2, 
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b) 

FIG. 2: Diagrams for p° — + 7r°7 decay. Here "•" denotes all potential meson loop correction, bold 
solid lines and dash lines denote type I complete propagators of p° and uj, which is defined in fig. I, 
"• • •" denotes all potential high order diagrams. In addition, the bold bold solid line in fig. 2-b) is 
defined as complete propagators (not type I) of p° in fig. 2-a). Thus the chain apporximation in 
fig. a) corresponds to renormalization of mass and wave function of external line p°, and the chain 
apporximation in fig. b) corresponds to one of internal line uj. 

every "•" denotes contribution from all potential meson loops. Thus the renormalization of mass and wave function 
of p° and uj is necessary. It should be pointed out that every bold solid line or bold dash line , in fig. 2 is defined 
in fig. 1. It means that the effect of p° — uj mixing has been included completely via summing all diagrams of fig. 2. 

There are two kinds of contributions to p° — ► 7r°7 decay. They are shown in the fig. 2-a) and 2-b), and respectively 
correspond to nonresonant contribution and contribution of resonance exchange. Renormalization of mass and wave 
function of external line p° is present in both of fig. 2-a) and 2-b), while one of internal line uj is present in fig. 2-b) 
only. For external boson line, Dyson has shwon at adiabatic limit p2|, 

(q 2 - m 2 )A F (q 2 )\ q 2 =m i — » iz\ , (10) 

where A F (q 2 ) denotes chain approximation to the exact propagator, in which the mass renormalization has been 
performed. This relation hold all of external p , ir° and photon fields. Meanwhile, the renormalization of mass and 
wave function of internal line uj yields its complete propagator as follow 

a m M 2 ) = ?=• ( n ) 

q 2 - m 2 + iTu{q 2 )\/q 2 

Here the dynamical width T UJ (q 2 ) is generated by pion loops. It can be also determined by unitarity of S-matrix. 
According to the above discussion, the LSZ reducation formula (||) becomes 

< ir (k) 7 ( qi ),out\p (q 2 ),in >= if% 0l {q 2 2 )\ ql=m 2f va ?qi^lq 2a e^ (12) 

where the superscript "(c)" denotes non-perturbative coupling, 

Ac) ,*\ t ( i\ , n p"(9 2 )/iJo 7 (<? 2 ) 

q z - m£ + i V9 1 u\Q ) 

In eq. (p^|), j p o 7! a 1 {q 2 ) and /u W o 7 (g 2 ) denote renormalized form factors (momentum-dependent coupling), which include 
meson loop correction. Labeling them by a superscript "(0)" and taking q 2 to mass shell of p, eq. ( |l3| ) is just the first 
equation of eq. (||) . The similar discussion holds for the second equation of eq. (||) . So that we finish the prove for 
eq. (||) in effective field theory formalism. 

B. Mixed propagator approach 

Alternatively, there is another well-known quantum mechanics method to deal with the p — uj mixing problem: the 
approach of mixed propagator p^, p"4| . This approach was developed even before discovery of QCD. The vector meson 
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propagator is given by (Renard representation) 

D^{q 2 ) = j ^d 4 ^-^(0|T{V;(x)K(0)}|0) 



= D(9 2 )v+^P(0)-O(s))w» (14) 
where s = q 2 , and the propagator function -D(s) is written in the following way 

D < s > - T^WVy < 15 > 

For multi- vector-meson channels, W{s) is complex mass-square matrix with non-zero off-diagonal elements in general. 

In order to define physical states measured by experiment, let us consider the single vector meson resonance channel 
case first. The D^ v of eq. (|l4|) is now the ordinary propagator of a vector meson. The reaction amplitude for a process 
by the medium of this vector meson resonance reads 

M ~ J?D^J% = (J, ■ J 2 ) s _ 1 w{s) (16) 
where and J% represent some currents and q^J^ 2 = 0- The reaction probability is 

The meson resonance mass measured in the experiment is real and is determined by the location of the maximum of 
<r(s)— peak in real s— axis. Then the resonance mass is determined by the following equation 

^-\ S -W( S )\ 2 = (18) 

Consequently, the mass-square of the resonance, M 2 , is determined by the solution of above equation, i.e., 

1 d d 

s = [ReW - 1^ s ( w * w + WW*)](1 - -r^ReWy 1 = M 2 (s). (19) 

Now return to p a — oj two channel case, W(s) and hence M 2 (s) are complex and real 2x2 matrices respectively. 
The mass determination equation for physical resonance states should read 

det[s-M 2 (s)} =0. (20) 

The physical states then are the eigenvectors of the real mass-matrix M 2 . Following refs|H| Q and using Breit-Wigner 
approximation, we have 



\IV puJ {s) ml I -i^/sT UI J 

where m pi and m u , are masses of p and u> in isospin basis. Considering Imn pw is small and hence ignorablejl6|], we 
get 

ReW = ( < H ^ (S) ) 
\U pul (s) mi, J 

Then, in dReW/ds, only off-diagonal elements of dli pul /ds left. Generally, for a broad class of models H pul (s) at 
{p— , u— )resonance energy region can be determined by taking VMD-type p — u mixing Lagrangian L pu} = f pw p^ u ^^v 
( yv v — d^V v — d u V^, V — p,oS). It leads to U pu (s) — f pu s which satisfies U(s = 0) puJ = required by generic 
consideration in ref.[[l7|. Thus, we have 

d „ . n 



y/sT Pl U puJ (s) 



ds 



n pw = f puJ = -^U m 2 ~ 6.7 x l(T d « 1 =► 1 - (dReW/ds) ~ 1. (21) 
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where |n pw | ~ AOOOMeV Wn has been used in the estimation. Furthermore, noting T^/m^ <C 1, we have 



\-r-T-\fO.O\Fi-r-r O I ^ 



2 a(sTj) 



'/'/ 



where the off-diagonal elements of M 2 matrix represent the p — lo mixing in the isospin basis. The physical p and w 
are eigenstates of M 2 . In other words, M 2 matrix can be diagonalized by the unitary 2x2 matrix C: 

CM 2 C^=( m " \\ 
\ ml) 

where 

C -[v i J' ■ (23) 

Consequently, the solutions physical state condition of eq. ( po|) are follows 

\P° P ) = |p?>-«, (P° P \ = \P° P ) 1 

\u p ) = |wj)+J?|p?), H| = K)t. (24) 

Under this transformation, we have 

CWC^ = (25) 

where z p = m 2 — im p T p , z w = — im w r w , and T = n pa; — 77(2:^ — z p ) are all defined in the physical state basis. The 
propagator function in the physical basis D p reads 

d p (s) = cis-wy 1 ^ = {s-cwc ] y 1 

= ( {s-zp)- 1 (s-z p )-iT( S -z^\ 
[(s-z^Tis-Zp)- 1 (s-z^)- 1 J 

( D p D P TD P \ 

— 1 PP PP LOU) | 

~ \ D p TD P D p ) ' 



(26) 



For V — X-vertex (V = p, lo and X represents other particles), f p x denotes the corresponding form- factor, f p x 
and fy X denote the coupling constants in the physical basis and in the isospin basis respectively. Since generally 



(a) . 




FIG. 3: The relation between the form factor of V — X vertex and the corresponding coupling 
constants for the vertices. The black dots in the vertex denotes the form factor. The single thick 
(dash) lines denote p— propagator Dp P (uj— propagator D^ u ), and the thick-cross-dash lines (or 
dash-cross-thick) line denote the mixed propagator Dp^ (or D F p ). The thin lines are external lines 
of X-particles. 

D p u = D p p ^ 0, the form factor is different from the corresponding coupling constant, i.e., f p x ^ Jvx- Frarri Fig.l, 
we have 

Dvvfvx = Dy p f pX + D p u f p x , (27) 
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with 



fpx = /S " tix = & + ■ (28) 



IpX — JpX 'UojX JujX — JuiX ~ i J fry 

In terms of eqs. (|26| ) and (^8|) the time-like EM pion form-factor is given, in the p — ui interference region, by 
F (s) - 1 4- \f p D p f F + f p D p f F 1 

1 TT\OJ ' U pry ppJ PTT7T J 017 LOLOJ WITKl 



P7 pp-> pTTTT 1 JW7 
fP fP fP fP I fP fP , fP fP \rp 

^ . JpyJpTnr ^ Jlj^Juj-k-k ^ \J P7 J uj-ktt ' J ujjJ pTTir) 

S-Zp s-z u (s - z p )(s - z u ) 

(fPfP + fP fP )t i 

\J p'jj UJ7T7T 1 J ijj'y J ottttJ \ / 



1 + ifULr + P7 )(— (29) 



with 



£e** - [-7?- w ][1 + W _ 3/ ] , (30) 

O Zp Z^j Zp Zq; 

where / P7 = 3/i, 7 and f^n-n = have been used, and <j) is Orsay phase. Using W pu} ~ — 4000MeU 2 [|l7| in eq. (|3C]), 



we obtain that £ ~ 0.012 and <fi is equal to about 100° 101° as s varies from m 2 p tp m 2 . These predictions are in 

good agreement with experimental datajl5|, and hence the mixed propagator approach is legitimate to describe 
the p — ui mixing effects in the pion EM-form factor. 

Now, we study the anomalous-like p° — > ir°j and ui —>■ ir°j decays in terms of the mixed propagator approach. 

J£ AO) 

,2 I nm F ■' OT °7' 



. @, we 


have 




/ p o 7r o 7 — 


AO) 

J p 7T°7 


-77f (0) 


f F 







T f (0) 



m 2 - to 2 + im p T p 



/ p7r ° 7 ' (31) 



where eqs. (|24|) and ( |2 6|) have been used. Considering ~ 0.006 << 1 and T ~ n pw , we finally obtain eq. 
mixed propagator formalism. 



III. NUMERIC RESULT AND CONCLUSION 



The w — > 7r + 7r decay suggests that the on-shell amplitude IL^m^) is around -4000MeV 2 which is indeed very 



2 - to 2 + im w r w 



small for the most isospin conservation processes. However, for p° — > 7r°7 we can see that m 
(—18624 + 6577i)MeV 2 is also small due to narrow width of u>. Taking H pul (m 2 ) ~ — 4000MeV 2 and combining with 
eq. (||), we have / p o w o 7 //jo^ ~ 1.6 at large N c limit and chiral limit. Therefore, the hidden isospin breaking process 

p° — * ui — ► 7r°7 indeed plays significant role in p° — > 7r°7 decay. In addition, f U n°'y/f}o — 0.99 is obtained, so that 
the contribution from p° exchange in u decay can be omitted due to T p ^> and eq. (||). 

For predicting branching ratio for p° decay precisely, the precise on-shell p° — u> mixing amplitude is needed. The 
investigation of p° — lo mixing has been an active subject jy], |l3[ [l5| [iq| . In ref. (T^] the on-shell mixing amplitude has 
been determined as II pu (m 2 ) = — [(3500 ± 300) + (300 ± 300)i]MeVMn the pion form factor in the timelike region, 
but in this fit effect of "direct" unrir coupling is omitted. We quote this result as R\. Meanwhile, the ref. [jll| provided 
a complete theoretical study on this mixing, which included the effect of "direct" uTtit coupling, and is up to the next 
to leading order of N^ 1 expansion. The result is n pw (m 2 ) = -[(3956 ± 280) + (1697 ± 130)i]MeV 2 , and we quote 
this result as i?2- Then using B(p ± — > 7r ± 7) = (4.5 ±0.5) x 10~ 4 and B(u> — ► 7r°7) = (8.5 ±0.5) x 10~ 2 and assuming 

/ p o^o 7 = / P ±7r± 7 , we obtain @ 



^(p°-A)= 1 _ z„ (32) 



(11.05 ± 1.84) x 10" 4 for R x 

(12.25 ± 1.52) x 10" 4 for R 2 

Here we have used II pw (m 2 ) ~ m 2 II pw (m 2 )/m 2 . The above results strongly support the fit of the solution A in table 
1 instead of the solution B. In fact, from viewpoint of experiment, the solution A is also better than the solution B 
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(because of smaller \ 2 an d reasonable phase difference). Therefore, we conclude that the branching ratio for p° — > 7r°7 
is 

B(p° -► vr°7) = (11.67 ± 2.00) x 10~ 4 . (33) 

Here our phenomenological estimates (eq. (|32|)) are kept as references instead of averaging them and experimental fit. 

Finally, it is interesting to estimate the contribution of u> exchange in p° — > 777 decay. This process itself is isospin 
breaking due to electromagnetic interaction. The exact isospin symmetry implies f^L/fiSni = 3- Hence similar 
argument gives fpO^/f o ~ 1.06 and B(p° — > t]j)/B(lu — > 777) ~ 0.5. This result agrees with current experimental 
fits, B(p° -► 777) = (2.4+°;*) x 10- 4 and B{lu -► 777) = (6.5 ± 1.0) x 10- 4 [§. 

To conclude, we show that a hidden effect of isospin symmetry breaking plays essential role in p° — > ir°j decay. The 
transition amplitude is derived by effective field theory approach and mixed propagator approach respectively. The 
result yielded by two independent method is matched each other. Our result is also supported by some experimental 
evidences. It indicates that the current datum for this process in PDG should be corrected. 
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